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Abstract 

We show how Stokes' Theorem, in the fashion of the Generalised Cauchy Formula, can be applied for 
computing double-cut integrals of one-loop amplitudes analytically. It implies the evaluation of phase- 
space integrals of rational functions in two complex-conjugated variables, which are simply computed by 
an indefinite integration in a single variable, followed by Cauchy's Residue integration in the conjugated 
one. The method is suitable for the cut-construction of the coefficients of 2-point functions entering the 
decomposition of one-loop amplitudes in terms of scalar master integrals. 



Unitarity and analyticity are well-known prop- 
erties of scattering amplitudes pQ. Analyticity 
grants that amplitudes arc determined by their own 
singularity-structure, while unitarity grants that the 
residues at the singular points factorize into products 
of simpler amplitudes. Unitarity and analyticity be- 
come tools for the quantitative determinaton of one- 
loop amplitudes [2] when merged with the existence 
of an underlying representation of amplitudes as a 
combination of basic scalar one-loop functions [3J. 
These functions, known as Master Integrals (Mi's), 
are n-point one-loop integrals, I n (1 < n < 4), with 
trivial numerator, equal to 1, characterised by ex- 
ternal momenta and internal masses present in the 
denominator. Important improvements of unitarity- 
based numerical algorithms also make use of the gen- 
eral structure of one-loop integrands [H [5J [7] . In 
the context of unitarity-based algorithms, the issue of 
computing one-loop amplitudes can be addressed in 
two stages: the computation of the coefficients; and 
the actual evaluation of the Mi's themeselves. The 
principle of a unitarity-based method is the extrac- 
tion of the coefficients multiplying each MI by match- 
ing the multiparticle cuts of the amplitude onto the 
corresponding cuts of the Mi's. 

Cutting a propagating particle in an amplitude 
amounts to applying the on-shell condition and re- 
placing its Feynman propagator by the corresponding 
<5-function, (p 2 — m 2 + i0) _1 — ► (2ni) S^(p 2 — to 2 ). 
As a result, the original function is substituted by a 
simpler one, easier to compute, which, nevertheless, 
still carries non-trivial information. In fact, the n- 
particle cut of /„ appears in the O-trascendentality 
term (rational or irrational) of the corresponding 
cut-amplitude, multiplied by the same coefficient of 
I n in the decomposition of the complete amplitude. 
Higher-transcendent ality terms, such as logarithms, 
are associated to the cuts of higher-point Mi's. 

In general, the fulfillment of multiple-cut condi- 
tions requires loop momenta with complex compo- 
nents. Since the loop momentum has four com- 
ponents, the effect of the cut-conditions is to fix 
some of them according to the number of the cuts. 
Any quadruple-cut [5] fixes the loop-momentum com- 
plctly, yielding the algebraic: determination of the 
coefficients of I n ,(n > 4); the coefficient of 3- 



point functions, ^3, are extracted from triple-cut 
[3 EH El H2 03] ; the evaluation of double-cut 
[HlIliaQIlliniQIlini^ necessary for ex- 

tracting the coefficient of 2-point functions, Ii\ and 
finally, in processes involving massive particles, the 
coefficients of 1-point functions, Ix, are detected by 
single-cut [T2l I2TH |2"T] . In cases where fewer than 
four denominators are cut, the loop momentum is not 
frozen: the free-components are left over as phase- 
space integration variables. 

In this letter, we show a novel efficient method for 
the analytic evaluation of the coefficients of one-loop 
2-point functions via double-cuts. Spun-off from the 
spinor-integration technique [14 } IT5 l [TB I IP71 li"8l [19] . 
the method hereby presented is an application of 
Stokes' Theorem. We analyze the double-cut of mass- 
less particles in four-dimensions, which also is the es- 
sential ingredient for the phase-space integration in 
the general case of one-loop massive amplitudes in 
dimensional-regularization [TH1 [T71 [T5J Q3B IH] • 
Due to a special decomposition of the loop- 
momentum, the double-cut phase-space integral is 
written as parametric integration of rational function 
in two complex-conjugated variables. By applying 
Stokes' Theorem, the integration is carried on in two 
simple steps: an indefinite integration in one vari- 
able, followed by Cauchy's Residue Theorem in the 
conjugated one. 

The coefficients of the 2-point scalar functions, be- 
ing proportional to the rational term of the double- 
cut, can be directly extracted from the indefinite in- 
tegration by Hermitc Polynomial Reduction. 

In a framework where factorization properties of 
scattering amplitudes are accessed via complex mo- 
menta, the double-cut integration presented here can 
be considered as the natural extension of the tech- 
nique used to prove BCFW-recurrence relation for 
tree-level amplitudes [55]. In the latter case, scatter- 
ing tree-amplitudes are holomorphic functions, de- 
pending only on one complex variable, and Cauchy's 
Residue Theorem is sufficient for their complete de- 
termination. In the case of the double-cut of one- 
loop amplitudes, where the integrand depends on two 
complex-cojugated variables, Stokes' Theorem in the 
fashion of Generalised Cauchy Formula, becomes the 
driving principle. 
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1. Double-Cut 

- Phase-Space Parametrization. The starting point 
of our derivation is the spinorial parametrization of 
the Lorentz invariant phase-space (LIPS) in the K 2 - 
channel [2J El [15] , 



d 4 h S M (£i) S^dh-K) 2 ) 



K 2 



mil] r k 



obtained by rescaling the original loop- variable £^ as, 



(2) 



with £\ — £ 2 — 0. In terms of spinor variables, the 
rescaling reads, 



|4) = Vi \£) , \h] = Vt\l], 



(3) 



where t, the rescaling parameter, is frozen as a con- 
sequence of the (second of the) on-shcll conditions, 
and £^ becomes the new loop integration variable. 
- Change of Variables. We take two massless mo- 
menta, say p^ and fulfilling the conditions, 



< 2 = 0, 



P 



2p ■ q = 2p ■ K = 2q ■ K = K 2 , (4) 



and decompose £ ll in a basis of four massless momenta 
constructed out of them, 

z z 
t li =p ll + ZZ q^ + - (q\^\p\ + - (p\l»\q\ ■ (5) 

Notice that the vectors M^M an( j (phM are triv _ 
ially orthogonal to both p^ and q^. The above decom- 
position can be realized starting from the definition 
of 1^ in terms of spinor variables, £^ — ^ 7 2 ^ , and 
performing the following spinor decomposition, 



\£) = \p) + z\q) , \£] ee \p] + z\q] 



(6) 



By changing variables (\£), \£}) — > (z, z) as in ([6]), and 
using Q, one can write, 



(£ d£) [d£ £} = K 2 dz dz , 
(£\K\£] = K 2 (1 + zz) 



(7) 
(8) 




£i-K 

Figure 1: Double-cut of one-loop amplitude in the if 2 -channcl. 

Hence, the LIPS in reduces to the novel form, 

1 



dz dz dtt z 5[t 



(1+zz) 



, (9) 



where t is a positive quantity as assured by the argu- 
ment of the 5-function. 

- Double-Cut Integration. The double-cut of a 
generic n-point amplitude in the ^-channel is de- 
fined as 



A ee / d 4 $ A t [ cc ( 



At e (\ti),\ii]) , (10) 



where A^ ^ are the tree-level amplitudes sitting at 



the two sides of the cut, see Fig 
as in ([3]), and using expression ( 
has, 



After rescaling £\ 
for the LIPS, one 



AT c (t,\£),\£}) At%tM),\g\) 
dz I dz I t 2 dt 5[ t ' 



(1 + zz), 

t aL+aR At ee W)M AT e {\i)M , (ii) 

whera <xl,r parametrizes the scaling behaviour of 
A^r. The i-integration can be performed trivially, 
because of the presence of the (5-function. Then, by 
using the decomposition ([5j [6]), the double-cut be- 
comes a double-integral, 



(12) 



A = <p dz dz f(z, z) 



where / is a rational function of z and z. As such, it 
can be expressed as a ratio of two polynomials, say 
P and Q, 



f(z,z) 



A t [ ee {z, z) A% ee {z, z) P{z, z) 



(1 + zz) a ^+ Q s+i 



Q(z,z) 



(13) 
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with the following relations between their degrees, 
deg z Q = deg z P + 2 , deg z Q = deg- z P + 2 . (14) 
We remark that the double integration in z- and z- 



variables appearing in Eq.( 12 1 will be properly justi 



fied in Sec|2] For the moment, with abuse of notation, 
we simlpy denote it as a convolution of an indefinite 
z-integral and a contour z-integral, which are the ac- 
tual operations we are going to carry out. 
To begin with the integration, we find a primitive of / 
with respect to z, say F, by keeping z as independent 
variable, 



F(z,z) -- 
so that A becomes, 

A = * dz F(z,z) = 



dz f(z,z) 



dz F s 



(15) 



(16) 



where F s is a short-hand notation for dF/dz. Be- 
fore proceeding with the final integration on the z- 
variable, let us analyse the structure of F. Since F is 
the primitive of a rational function, its general form 
can only contain two types of terms: a rational term 
and a logarithimc one, 



F(z,z) = F lat (z,z)+F lo s(z,z) 



(17) 



It is important to notice that the presence of the term 
F rat depends on the powers of t in Eq. (11): F rat 
can be generated, after integrating / in z, only if 
olr + oll > 0. The z-integration will be performed by 
applying Cauchy's Residue Theorem, therefore the 
final structure of the double-cut is determined by the 
nature of F. Namely, the z-integration of F rat [F log ] 
is responsible of the rational [logarithmic] term of A. 
- 2-point Function. We also know apriori that 
the double-cut of a 2-point scalar function in 4- 
dimension is a rational, or better (to account for the 
massive case as well) a non- logarithmic term; while 
the double-cut of higher-point scalar functions might 
contain logarithms (with if 2 -dependent argument). 
Hence, the coefficient of a 2-point function in the K 2 - 
channel will appear in A rat , the integration of F rat 



dz F rat (z,z) 



(18) 



where the z-integration is performed via Cauchy's 
Residue Theorem. The integrand F rat is rational in 
z, and contains poles whose location in the complex 
plane is a unique signature of the Feynman integral 
they come from [T71 [TH1 [21] ■ The choice of p and q 
specified in Eqs.Q grants that there exists a pole at 
z = associated to the 2-point function in the K 2 - 
channel, l2(K 2 ); while the reduction of higher-point 
functions that have l2(K 2 ) as subdiagram can gener- 
ate poles at finite z-values. Because of the presence 
of z, through the term (1 + zz), F lat is non-analytic. 
The Residue Theorem has to be applied by reading 
the residues in z, and substituting the corresponding 
complex-conjugate values where z appears. There- 
fore, the result of A rat can be implicitly written as, 



A rat = 27nS(Res 2=0 F rat (z,z) +Res z ^ F 1 ' at (z,z) 



(19) 



- Double-cut of the Scalar Function 7 2 - Let us eval- 
uate the double-cut of the 2-point scalar function I 2 , 
which also is a prototype example: 



Ah 



d 4 ^ 



dz dz it 2 dt 51 t 



(1 + zz) 



dz I dz 



(1 + zz) 2 



= <b dz 



(-1) 
(1 + zz) z 



(20) 



where we notice that the primitive in z, called F in 



(15 1, has only a rational term (the logarithmic con- 



tribution is absent), 

F(z,z) = [dz f(z,z) = = F rat (z,z). (21) 

J (1 + zz) z 

For the last integration in z, by applying the Residue 
Theorem, we take the residue of the unique simple 
pole at z = 0, since the term (1 + zz) = (1 + |z| 2 ) 
never vanishes, being always positive. The final result 
of the double-cut of the scalar 2-point function reads, 

AI 2 = (2m) Res z=0 ^ rat (z,z) = -2ni . (22) 
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- Coefficient of the 2-point Function. The expres- 
sion of the 2-point coefficient can be finally obtained 
by taking the ratio of A rat in ( 19 ) and the double-cut 



of I 2 in (22), 



C2 



A7 2 

-Res z=0 F rat (z,z) 



Res z#0 ^ rat (^) • (23) 



- Hermite Polynomial Reduction. To optimize the 
integration algorithm, one can use the so called Her- 
mite Polynomial Reduction (HPR), a technique en- 
abling the direct extraction of the rational term of 
the primitive of a rational function, without comput- 
ing the integral as a whole. Based on the square-free 
factorization of the integrand, HPR can be used to 
write the result of any integral of a rational function 
as a pure rational term plus another integral that, if 
explicitly computed, would generate the logarithmic 
remainder. 



As written in Eq. ( 23 ) , the coefficient of the 2-point 
function comes only from the term F rat , and not from 



F log , see Eqs.( 15 17 1; F rat is the rational term in the 



result of the z-intcgration of /, which is rational in z, 
see Eq.( 13 ). Therefore HPR is suitable for extracting 
_F rat out of the z-integration. 

The integration algorithm of Sec [I] can be imple- 
mented with S@M [25 together with the routine [25] 
for Hermite Polynomial Reduction. 



2. Stokes' Theorem 

In this section we give a formal definition of the z-z 
integration used in Sec|T] as an application of Stokes' 
Theorem for differential forms. In what follows, we 
use the notation: g z = dg/dz and g g = dg/dz. 

Let us recall that the complex 1-form 



1 



X 



-dz 



(24) 



z- z 

which is defined for all z except zq 7 is a closed form, 
1 \ ■ (-1) 



z ■ 



-0 



A dz 



(z - z y 



: dzAdz = . (25) 



We consider any complex smooth function T and dif- 
ferentiate the 1-form lu = J-\, 



(z — zq) 1 J-dz , 



(26) 



obtaining the 2-form, 

duj = dT A\= {z - zq)- 1 ^ dz Adz . (27) 

Now we take a domain D in the complex plane and 
apply Stokes' Theorem to doj. Due to the singularity 
of w at zo, we remove a tiny disk D(zq; r), centered at 
zq with radius r, from D. Then oj has no singularity 
in the regulated domain D r = D — D(zq; r), and we 
may apply Stokes' Theorem: 



dto 



dD r 



()D 



dD(z ;r) 



(28) 



Here dD(zo; r) is a circle 7 around the point zq, which 
is described by the parametric equation ^(t) = z + 



Since T{z^ 



') converges to T{zq) as the 



radius r shrinks to 0, the last integral in Eq.^28 



/ w = * / 

JdD(z ;r) JO 



F{z + re lt )dt 



(29) 



converges to 2nij r (zQ) as r goes to 0. Letting r 



in Eq.(28l, the disk D(zo;r) disappears and D r 
fills up D. Consequently Stokes' Theorem can be 
reformulated as, 



dui — / u> — l-niTizo) 

JdD 



(30) 



By using the explicit expression of u> and duj, in 



Eqs.(26 27 1, and rearranging terms, we obtain the 



so called Generalised Cauchy Formula or Cauchy- 
Pompeiu Formula, 



2-kiT(zq) 



Hz) 



dz 



dD z — Zq 



J z 



D z Zq 



-dzAdz. (31) 



Let us discuss two special cases. 

First, when T is analytic, T% — 0, hence we obtain, 



f(zo) 



1 

2~7ri 



Hz) 

dD Z — Zq 



dz 



(32) 
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which is the well-known Cauchy Formula, where dD 
is any closed curve surrounding Zq. 
Secondly, when T vanishes on the boundary of D, 
that is T\qd 



0, Eq.(31 1 becomes, 

— [[ Tl 

2ni 



dz A dz. 



D 



ZQ 



(33) 



where we used dz A dz = —dz A dz. 



The expression (|33|) is what needed to define prop- 

which 



erly the double-cut A given in Eqs.(12 
we rewrite here as, 
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A = I J f(z, z) dz Adz 



-dz Adz , (34) 

Id z ~ z q 

by identifying / = F s = F s /(z-z ), and F = 
J-/(z — zq), where the functions / and F were de- 
fined in Eqs. ( fl3| [l5{ |l7| ) . The integration domain, D, 
is the whole complex plane. The vanishing of J- on 
the boundary is granted by the structure of the ratio- 



nal integrand and relations ( 14 1 among the degrees 



of numerator and denominator. 

To deal with the general case, where more than one 
pole might appear, the calcualtion of A trivially gen- 
eralises, by the superimposition principle, to the sum 
of the residues at all the poles in z, 



A EE 



Fz dz A dz = 



-dz A dz 



E 



2m J2 ^fe), 

j'Gpoles 



(35) 



due to the subtraction of a disk around each of the 
z-poles from the domain D. 



Finally, Eq.(35l validates Eq.(19), hence the expres- 
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